Absence of mobility edge for the Anderson random potential on 
tree graphs at weak disorder 



Michael Aizenman 1 and Simone Warzel 2 

1 Departments of Physics and Mathematics, Princeton University, Princeton NJ 08544, USA 

2 Zentrum Mathematik, TU Miinchen, Boltzmannstrasse 3, 85747 Garching, Germany 



o 

(N 
Oh 

C/3 



PACS 73.20.Jc 
PACS 73.20.Fz 



Derealization processes 
Weak or Anderson localization 



Abstract -Our recently established criterion for the formation of extended states on tree graphs 
in the presence of disorder is shown to have the surprising implication that for bounded random 
potentials, as in the Anderson model, there is no transition to a spectral regime of Anderson 
localization, in the form usually envisioned, unless the disorder is strong enough. 
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Introduction. — Since the early studies of Anderson 
c~*| localization, random operators on tree graphs have pro- 
"£^vided a testing ground for insights on the effect of disor- 

Sder on quantum spectra and dynamics [l}j3j. The basic 
phenomenon to which this refers is the transition in the 

spectra of Schrodinger-type operators with random poten- 
tial from regimes of energies with only localized states, to 
energies with extended states which facilitate conduction, 
e.g. in the sense depicted in Figure [l] 

(N Our purpose here is to report on rigorous results which 
challenge some of the generally shared picture of the mo- 
0*N bility edge in such systems on the Bethe lattice, i.e. within 
C*a regular tree graph (whose degree here is K + 1). These 
results are not expected to affect the picture of Ander- 
• .son localization for one-particle Hamiltonians in finite 
^dimensions. However, in view of analogies which were 
drawn between tree graphs and many-particle configura- 
!_i tion spaces [4j (cf . [5j|6] ) the mechanism which causes the 
^surprising effect may be worth paying attention to. Its 
root cause is the formation of extended states through 
fluctuation-enabled resonances between states which ini- 
tially (i.e., up to a certain scale) may appear to be local- 
ized. 

As was explained in [7], the exponential increase in the 
volume of the relevant configuration space implies that 
resonances whose likelihood decays exponentially in the 
distance will nevertheless occur, and predictably so, pro- 
vided the rate is small enough. Somewhat heuristically: 
states which may locally appear to be localized have arbi- 
trarily close energy gaps (AE 1 ) with other states to which 
the tunneling amplitudes decay exponentially in the dis- 




Fig. 1: A model setup for quantum conduction through the 
graph (after [24]): particles are sent at energy E = k 2 + [/wire 
down a wire which is attached to the graph at x = 0. The 
reflection coefficient is strictly less than one, |-Re| 2 < 1, exactly 
if (JsJ> holds, which is also the condition for the presence of an 
absolutely continuous component in the spectrum of H\(uj). 



tance R, as e~ Lx ^ E ' R . Mixing between two levels would 
occur if AE <C e~ Lx ^ R . The latter is a rather stringent 
condition, however since the volume grows exponentially 
fast (as K R ) extended states may form. A sufficient con- 
dition for that is: L\(E) < log if (for the exact condition 
one needs to take into account also the effects of large 
deviations) . 

As a consequence, even under circumstances which lo- 
cally seem to imply localization, e.g. as in the Lifshitz tail 
regime of low density of states, we find extended states. 
Some basic rigorous results related to this mechanism were 
presented the recent work |8j. Our focus here is on the 
somewhat surprising implication for bounded random po- 
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tentials, as featured in a class of basic models of the An- 
derson transition. 

The basic model. — A specific example of a one par- 
ticle Hamiltonian to which our discussion relates is 



1 a * a y 



xeg 



(1) 



where the first sum is over pairs of neighboring sites of a 
regular tree graph Q, i.e., the Bethe lattice, and a\a are 
fermionic creation and annihilation operators. The dis- 
order is expressed in the random 'onsite energies' e x (uj) 
which are given by independent identically distributed 
(iid) random variables, e.g., of a uniform distribution in 
the interval [— W/2, W/2]. A more general class is de- 
scribed below. 

In the one-particle sector, the above Hamiltonian may 
also be written, as is often done in the mathematical lit- 
erature on the subject, as 



H\(u) = T + \V{lo) on£ 2 (g), 



(2) 



where: T is the graph's adjacency operator, which acts as 
(Tip)(x) = Y^ y ~x i s a random potential, and u 

represents the randomness. The strength of the disorder is 
controlled by the parameter A (> 0), which in the notation 
of corresponds to W/(2t). 

The results presented below are formulated for random 
potentials whose values at the graph sites, x £ 0, form 
iid random variables V x (u>) whose probability distribu- 
tion is of the form p(V)dV with a density function p(V) 
which satisfies a simply stated regularity condition, which 
is spelled below. Our emphasis here is on bounded ran- 
dom potentials, and more specifically on the case when 
the support of p, which is the closure of the set on which 
p{V) 7^ 0, is exactly [—1, 1]. A simple example to which 
the discussion applies is the 'Anderson random potential' 
with p(V) constant over the interval [—1, 1] and zero else- 
where. 

On homogenous graphs, such as the regular lattice Z d or 
the homogenous trees discussed here, the spectrum itself 
is easy to determine for random operators such as H\{ui): 
by ergodicity arguments [9 -11 it is the set sum of the 



spectrum of T and of XV. More explicitly, for tree graphs 
with branching number K (both the rooted and the ho- 
mogeneous case) and potentials described above, the spec- 
trum of H\(uj) is for almost all realizations given by the 
non-random set 



a(H x ) = [-\E x \,\E x \], with£ A 



(2VK + A) (3) 



(An auxiliary result of [12] clarifies that the spectral mea- 
sures associated to the different vectors S x £ £ 2 (G) are 
almost surely all equivalent.) Thus, as the strength of the 
disorder is increased from A = upward the outer edges 
of the spectrum spread at the linear rate A. However, the 
nature of the spectrum is somewhat less obvious. 



Localization 
(p.p. spectrum) 




-2VK 2\/K 



Fig. 2: A schematic depiction of the previously expected phase 
diagram for the Anderson model on the Bethe lattice (the solid 
line) and the correction presented here (dashed line). The pre- 
vious incomplete theoretical analysis |2[[3| has received some 
support from extensive numerical work [25] . Our analysis sug- 
gests that at weak disorder there is no localization and the 
spectrum is purely ac. While the proof of that is incomplete, 



we do prove that for A < (y/K - 1)72) there is no mobility 
edge beyond which localization sets in. 



The phase diagram. — It is known that in one di- 
mension complete localization sets in even at arbitrarily 
weak disorder, with the localization length being finite 



(except at isolated energies in some special cases 13 14]) 



and dependent on the disorder strength 15 16 . However, 
it is expected 17 that above two dimensions at weak 
and intermediate disorder the random operator exhibits 
both spectral regimes, with a phase diagram, which in 
case supp p = [—1,1] is roughly outlined in Fig. [2[ 

The above picture has also been expected to describe 
the phase diagram for the Bethe lattice. That is: 

i. At weak and moderate disorder a mobility edge has 
been expected, beyond which the spectrum consists 
of a dense countable collection of eigenvalues with lo- 
calized eigenfunctions. In the central regime the spec- 
trum is absolutely continuous (a property henceforth 
referred to as ac) with extended states which play an 
important role for conduction. 

ii. The extended states disappear at strong enough dis- 
order (A > X s d(K)), where complete localization pre- 
vails. 

Our results show that for the Bethe lattice this picture 
needs to be modified. 

Significant parts of this picture are supported by rig- 
orous results, in particular at strong disorder [7, 18 20 
(though some questions remain as to the precise asymp- 
totics of X s d(K) for K — > oo). The existence of ac spec- 
trum for tree graphs was also established [7|[8)[2l}|23] . 

Numerically, both Q and the above point i. have 
been hard to see in exact diagonalization of finite sys- 
tems, and/or in simulations of the infinite Bethe lattice. 
Researchers who carried such work report that Eq. |3| 
(whose proof is undisputed, and not complicated) is hard 
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to reconcile with numerical evidence 24 25 . The reason 



is easy to understand: since at the spectral fringes the 
mean density of states is very low, large simulations are 
needed in order to find evidence of the existence of eigen- 
states there. Naturally, it should be even harder to probe 
reliably the nature of states with energies in that regime. 
These eigenstates are associated with rare occurrences of 
nearly extremal values of the potential, realized over suffi- 
ciently large connected sets. Since such events occur only 
in well separated locations, the states may initially ap- 
pear to be localized. It is perhaps because of that, as well 
as the general theoretical guidance offered by the earlier 
works on localization, that some numerical studies have 
been interpreted as providing evidence for the existence of 
a mobility edge which connects to the unperturbed spec- 
trum o~(T), as sketched by the solid line in Figure [2j 

The general picture of the phase diagram for unbounded 
potentials was only recently resolved. As described in 
Refs. [7||8], on tree graphs the ac spectrum is fairly ro- 
bust, and it even appears under arbitrarily weak disorder 
at energies which are remote from that of the unperturbed 
spectrum cr(T) = [— 2\K, 2\K\. The mechanism which 
was identified there is also behind the more surprising re- 
sults for bounded potentials to which we turn next. 

Our main result. — It is helpful to first introduce 
the following quantity, which (for reasons related to a dy- 
namical system perspective) is referred to as the Lyapunov 
exponent: 



L X {E) 



limE 



log 



1 



'H x -(E + i V ) 



|0> 



, (4) 



where the expectation value sign E (•) stands for average 
over the disorder. As is explained in greater detail in j8j, 
on tree graphs the Lyapunov exponent provides the typical 
decay rate of the Green function: 



l H x -(E + iO) 



Const, e 



-Lx{E)\x\ 



(5) 



with (E + iO) indicating the limit r/ \ 0, as in Q, which 
exists for almost all energies. 

The results mentioned below are derived under a 
regularity conditions for the probability distribution 
of the potential. In its general form (a simpler one is 
presented below) it is stated in terms of the "minimal 
function" , which for an integrable function g on R, we 
define as: M g (V) := inf„ e (o,i] ^ u _ v ^ u g{u)du/(2v). 

Definition We call probability distribution on the line Ir- 
regular if it is absolutely continuous and its density func- 
tion is bounded relative to its minimal function, i.e., 



p(V) < bM p (V) , for almost every V € 



(6) 



at some b < 00. 

A simple sufficiency condition for (|6| is that p is bounded 
and there is a partition of the line into a finite collection of 



intervals within each of which p is continuous and mono- 
tone, except possibly at the interval's boundary (c.f. [8]). 

The results presented here start from the following ex- 
tension of the main finding which was reported in [7] , and 
which was now shown to be valid also for bounded random 
potentials [8]. 

Theorem 1. Let V be a random potential whose values 
are given by independent identically distributed random 
variables with an M-regular probability distribution. Then, 
for the operator defined by ^ on a K -regular tree, at al- 
most every energy E 6 o~(H x ) at which also 



L\(E) < log if . 



(7) 



one has, with probability one: 
1 



Im (x\ 



—r \x) > 0, at any x € Q . (8) 



1 H x - (E + »0) 

The dynamical significance of ^ was nicely illustrated 
in |24| : if quantum particles are sent coherently at energy 
E = k 2 + U W i re down a wire attached to the graph at the 
vertex x, the reflection coefficient Re satisfies \Re\ < 1 if 
and only if (|8) holds (cf. Fig. [I] and '26'). More explicitly, 
in the wire the particles' state is described by plane waves, 
the wave function at position £ being e lkE ^ + Re e b S 
The stationary state is given by a solution of the equation 
(H-E)^f = in £ 2 {g)®L 2 {Z + ) (with Q the graph and Z + 
the wire). The matching conditions for this equation at x 
relate the reflection coefficient Re to the Green function's 
value at (x,x), and yield the equivalence: 



\Re\ < 1 if and only if Q is satisfied. 



(9) 



The condition (|8| is also of direct spectral significance: 
the quantity 7r _1 Im (x| Hx _^ E +io) \ x ) 1S the density of the 
ac component of the spectral measure associated with the 
state \x) . Thus the criterion of Theorem [T] readily implies 
that in any energy interval within a(H x (u))) in which ([8]) 
holds on a set of energies of positive Lebesgue measure, the 
random operator almost surely has some ac spectrum (and 
hence also extended states, at the corresponding energies) . 

It is therefore of relevance to learn about the region in 
the phase diagram in which the condition ^ holds. In the 
absence of disorder, i.e. at A = the Lyapunov exponent 
is easy to compute, and one finds: 

= log Vk 
L (E) { G (log VK, log K 
> log if 



\E\ < 2VK, 
2VK < \E\<K- 
\E\ > K + 1. 



(10) 

Notably, the range of energies at which the condition |7| 
holds at A = (that is \E\ < K + I), is strictly larger 
than the spectrum of iio (which is \E\ < 2\K). In the 
context of unbounded potentials (for which the spectrum 
a(H x (uj)) covers the whole line as soon as A ^ 0), em- 
ploying continuity arguments we previously showed that 
Theorem [l] carries the following implication [t[[8] . 
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Corollary 1. Under the assumptions in Theorem^ for 
an unbounded random potential, with supp p = R, for any 
interval I C (— (K+l), K+l) at sufficiently small disorder 
the random operator has ac spectrum within I . 

For bounded random potentials the spectrum of H\(uj) 
expands at a bounded rate as A is increased from 0. In that 
situation the above continuity statement does not yield 
much information about the spectrum outside of o~(T). 
For the same reason, the general result of jT9] which in 
essence proves localization at weak disorder for energies 
at which (0\(H — E)^ 1 )^) is summable, is also not appli- 
cable. To apply Theorem 1 to this case, we present here 
the following result. 

Theorem 2. Let H\(lu) be a random operator of the 
form ([2| with a bounded random potential taking values 
in [—1, 1]. Then for all X small enough, and in particular: 



X<A K :=(VK-l) 2 /2, (11) 
the condition Q holds for all energies satisfying 

\E\ > \E X \ - 5(X) (12) 

at some (5(A) > 0. 



One may note that the regime defined by ( 12 ) includes 
the spectral edges 



E\, Ex 



6(X)) U (\E X \ - S(\),\E X \\ 



(13) 



as well the complement of the spectrum. The current proof 
yields only a rapidly vanishing value for (5(A). However, 
as mentioned in the discussion below, the result may be 
indicative of a much stronger statement. 

This leads us to the following somewhat surprising con- 
clusion. 

Corollary 2. Let H\(uj) be a random operator of the 
form ([2| with a bounded random potential whose density 
function p(V) satisfies the regularity assumptions of The- 
orem^ and has is supported in supp p = [—1,1] ■ Then, 
for A small enough, as in Theorem 2, near the spectrum 's 
boundary (i.e., in the spectral edges described by (13|J the 



random operator almost surely has only purely absolutely 
continuous spectrum. 

The proof is by a direct application of Theorems 1 and 
2, combined with the criterion of [27,28|, to which our at- 
tention was called by Mira Shamis. The latter states that 
for operators with random potential a sufficient condition 
for the spectral measure to be purely absolutely contin- 
uous in a given interval is that for almost every energy 
there condition ^ holds almost surely. 

Thus, for the bounded random potentials discussed 
here, at weak disorder there is no mobility edge beyond 
which the spectrum is pure point and localization sets in 
as has been expected (cf . Fig. [2]) . 



Outline of the proof of Theorem 2. First let 
us explain the explicit condition (111. Underlying it is 



the following deterministic lower bound, which is valid for 
any potential satisfying |U(x)| < 1 at all x, and E\ = 
-2-JK-X: 



1 



1 Ho + XV - E : 



|0> > <0| 



1 



H Q - (Eo - 2A) 



10) • (14) 



This relation (between positive terms) follows by the 
monotonicity of the operator- valued function A \— > 1/A, 
for positive operators A, combined with the elementary 
observation that: 

< Ho + XV -E x < H Q -{E Q - 2A) (15) 

Hence, for any A < [(K + 1) - 2y/K]/2 = A K : 

Lx(Ex) < L Q (E -2X) < log if. (16) 



where the second inequality is by ( 10 ). Consulting the ex- 



act values of the Lyapunov exponent, which is computable 
for A = (cf. M), one can say more: 

foranyA<A K : L x (Ex) <\ogK - C(A K - X) , (17) 

at some C > 0. 

Thus, for A < Ak the condition Q continues to hold 
at the boundary of the spectrum. To prove Theorem [2] 
we need to show that it continues to do so also within an 
interval which extends into the spectrum. Following is a 
summary of the main ingredients of the proof of that. 

1. To study the resolvent Gx(x,y;E) = (x\ Hx _ 1 E _ i0 \y) at 
energies in the range E G [E x , Ex + AE] , we compare it 
to the corresponding quantity for the restrictions of the 
operator to finite subtrees of depth R (measured from the 
root). We denote the corresponding operator by H x (ui), 

( Ft) 

and the matrix elements of its resolvent by G\ (x,y;E). 

2. It is of relevance here to note that by the variational 
principle the restriction raises the lower edge of the spec- 
trum. The values of AE > and R < oo are selected 
so that with probability sufficiently close to one (as deter- 
mined by the rest of the argument) 



inf o-(H[ R \u)) >E x + AE. 



(18) 



This is done with the help of the explicit probability esti- 
mate 

P (inia{H ( x R) ) <E X + AE^j < C K R {AEf' 2 (19) 

at some C < oo. This is based on a familiar argument for 
such a purpose, using as input the available bounds on 
Etre -tff (The bound seems far from optimal, as it does 
not take into full account the expected Lifshitz tail decay 
of the density of states at the edge by which the power law 
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(AE) 3 / 2 in ( 19 ) could be replaced by a much faster decay- 
ing rate, e.g., perhaps as discussed in the recent work [29]). 



configurations when cither (18) or (23) fails is controlled 



3. Proceeding under the assumption that (18) holds, 



we show that the finite volume analog of the Lyapunov 
exponent, L x R \e), for which H\(uS) is replaced in Q 

by H x R \u), satisfies the analog of ^ for all energies 
E < E\ + AE. The argument is facilitated by positivity 
arguments which are applicable under the above assump- 
tion. 

4. The difference \L X R) (E) - L\(E)\ is to be estimated 
using bounds on the difference of the corresponding Green 
functions. For the latter, the resolvent identity yields 



\g x (0,0; E)-G (R) (0,0; E) 



< 



\G{ R H0,x-;E)\\G x (0, X ;E)\ =: S{ r \e) . (20) 

x: \x\=R 

with \x\ the distance to the root, and x_ the site preceding 
x relative to the root. 

5. To estimate the error term s{ (E), we prove that 



in case (18) holds \G R (0, x_; E)\ decays at a faster rate 



than 1/Vif . More explicitly, for all A > there are 

(5(A) > 0, and C(A), C(\, s) < oo such that: 

i. Except for a rare event, of probability e(R, A, A) 



g{ R) (0,x.;E)\ < 



C(A) 



K (l/2+S(\))R 



(21) 



simultaneously for all x with |x| = R. 

ii. For the other term the decay of the moments s E (0, 1) 

is not slower than at the rate 1/VK , in the sense that 
for x as above 



E(\Gx(0,x;E)\ 
As a consequence, the event 



< 



C(X,s) 

J(s\x\/2 



(22) 



S{ R) {E) < Go(0,0;E -2\)K- 5 ^ R / 2 (23) 

with a suitably small <5(A) > has a probability which 
decays exponentially in R. 

6. The above results are combined into the observation 



that in case <jl8h and (23) hold: 



\G x (0,0;E)\ > G {R) (0,0;E) - S (r \e) (24) 
> G o (0, 0; E a - 2A) (1 - K R e 2RL ^ E «- 2 ^ - K' 5 ^ 2 ) . 

where the last factor can be made arbitrarily close to 
1 through an appropriate choice of R. The contribu- 
tion to the expectation of value of log |Ga(0, 0; E)\, a 
quantity which is of finite variance, due to the potential 



through the square root of the probability of that rare 
event. Since the main contribution is described by 
( [24| we find that L\(E) is bounded from above by 
L (E a — 2 A) = -logG (0,0;i?o — 2 A) plus a correction 
which vanishes as AE — > 0. This concludes the proof of 
Theorem [2 

The full details of the analysis outlined above are pro- 
vided in the lengthier and more technical article 8 . Our 
goal in the above summary has been to convey the essence 
of the result, and the idea of the proof. Since the focus 
here was limited to what is presently provable rather than 
what may be more generally true, we add comments on 
the latter in the discussion which follows. 

Discussion. — The relevant mechanism behind the 
results reported here is the formation of extended states 
out of local modes which up to a certain distance may ap- 
pear to yield localized eigenstates. These may join to form 
extended states through rare resonances, when the gap be- 
tween the local self energies is smaller than the tunneling 
amplitude. The latter decays exponentially in the dis- 
tance. Therefore the two important factors for the effect 
are the local fluctuations in the self energy, and the expo- 
nential growth of the volume of the configuration space. 
Under the Lyapunov exponent condition ^ the rare reso- 
nance conditions are actually met even at large distances. 
The implication for the bounded potential, namely the 
absence of a mobility edge at weak disorder at the edges 
of the spectrum which is reported here, is then a natu- 
ral consequence of the (at least) partial continuity of the 
Lyapunov exponent. 

Let us add that more may be true: Guided by the expec- 
tation that L\(E) is jointly continuous in (A, E) through- 
out the spectrum, we conjecture that for bounded poten- 
tials, of sufficiently regular distribution, at weak enough 
disorder there is no localization at all and the spectrum is 
purely absolutely continuous. Nevertheless, both the gen- 
eralized eigenstates and the time evolution may given the 
impression of the existence of localization, but that should 
be limited to finite time and distance scales . 

Another question which may warrant further attention, 
is whether the mechanism discussed here plays a similar 
role in case of particle systems with short range interac- 
tions. There are profound differences between the two 
systems, but also there the volume of the set of configura- 
tion reachable by R steps (when R is of the order of the 
volume) grows exponentially fast. 
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